WO 2005/058156 PCT/DK2004/000722 

1 

A system and a method for analysing ECG curvature for Long QT 
Syndrome and drug influence. 

The present invention relates to a system for analysing drug influence on ECG curva- 
ture and Long QT Syndrome where at least one among a number of different parame- 
ters is isolated, which system has an input means connected to an ECG source, where 
5 the different parameters of a received ECG curvature are indicated and/or isolated and 
for indicating possible symptoms. 

The present invention further relates to a method for analysing drug influence on ECG 
curvature, which curvature contains a number of parameters. 

10 

The heart generates an electrical signal which can be measured as an ECG, which can 
be recorded as an ECG diagram. The waves in the ECG-signal P,Q,R,S,T and U are 
due to depolarisation and repolarisation of the heart. 

1 5 Intervals and complexes are illustrated on a typical ECG curvature, see fig. 1 in order 
to illustrate the different curve sections isolated by the analysis and referred to by pa- 
rameters, intervals and complexes comprising the following annotations, where 

Ponset 2 marks the beginning of the P wave. 

20 

Qonset 4 marks the beginning of the Q wave. 
Rpeak 6 marks the top of the R wave. 
25 The Jpoint 8 marks the end of the S wave. 

Tstart 10 marks the beginning of the T wave. 
Tpeak marks 1 2 the top of the T wave. 

30 

Tend 14 marks the end of the T wave. 

The QT interval starts at Qonset 4 and ends at Tend 14. 

35 

The QT curvature is the part of the ECG curvature between Qonset 4 and Tend 14. 
The RR interval goes from one R-peak 6 to the following R-peak 7. 



CONFIRMATION COPY 



WO 2005/058156 PCT/DK2004/000722 

2 

US 5 9 749,367 describes a heart monitoring apparatus and method wherein an electro- 
cardiograph signal is obtained from a patient and processed to enhance the salient fea- 
tures and to suppress noise. A plurality n of values representative of the features of the 
electrocardiograph signal are generated and used in a Kohonen neural network to gen- 
5 erate an n dimensional vector. This vector is compared with a stored plurality m of n 
dimensional reference vectors defining an n dimensional Kohonen feature map to de- 
termine the proximity of the vector to the reference vectors. If it is determined by the 
Kohonen neural network that the vector is within or beyond a threshold range of the 
reference vectors, a signal is the output, which can be used to initiate an event such as 
10 the generation of an alarm or the storage of ECG data. 



US 2002/143263 describes a system comprised of a medical device and a method for 
analyzing physiological and health data and representing the most significant parame- 
ters at different levels of detail, which are understandable to a lay person and a medi- 

1 5 cal professional. Low, intermediate and high-resolution scales can exchange informa- 
tion between each other for improving the analyses; the scales can be defined accord- 
ing to the corresponding software and hardware resources. A low-resolution Scale I 
represents a small number of primary elements such as the intervals between the heart 
beats, duration of electrocardiographic PQ, QRS, and QT-intervals, amplitudes of P-, 

20 Q-, R-, S-, and T-waves. This real-time analysis is implemented in a portable device 
that requires minimum computational resources. The set of primary elements and their 
search criteria can be adjusted using intermediate or high-resolution levels. At the in- 
termediate-resolution Scale II, serial changes in each of the said elements can be de- 
termined using a mathematical decomposition into series of basis functions and their 

25 coefficients. This scale can be implemented using a specialized processor or a com- 
puter organizer. At the high-resolution Scale III, combined serial changes in all pri- 
mary elements can be determined to provide complete information about the dynamics 
of the signal. This scale can be implemented using a powerful processor, a network of 
computers or the Internet. The system can be used for personal or group self- 

30 evaluation, emergency or routine ECG analysis or continuous event, stress-test or bed- 
side monitoring. 
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The aim of the invention is to achieve a system and a method for diagnosing Long QT 
Syndrome in an objective, fast and effective way by indication of a number of symp- 
toms derivable from an ECG curve. A further aim of the invention is to achieve an 
effective test of drug influence on ECG curvature. 

5 

This can be achieved with the system previously described if a first number of selected 
parameters from at least three main groups, which groups comprise parameters of 
symmetry, flatness, duration and/or complexity, which parameters are combined in at 
least a first mathematical analysis, which relate to or are indications of certain dis- 

1 0 eases, where said diseases are known to influence the ECG curvature, where the result 
of the analysis can be represented as a point in a coordinate system comprising at least 
one axis where the system can compare the actual placement in the coordinate system 
with a number of reference parameters stored in the system for indicating symptoms 
or diseases having influence on the ECG curvature, where the system analyses the QT 

1 5 curvature of the ECG curvature for indicating Long QT syndrome. 

Hereby, it is achieved that any symptom of Long QT Syndrome having an indication 
(influence) in the ECG curvature can be detected in an objective, automated and very 
fast way. The system might be used under field conditions such as in ambulances or in 

20 other situations where a fast indication of heart diseases is needed in order to help the 
patient in a correct way as early as possible. The analysis that takes place in an ambu- 
lance on its way to the hospital can by transmitting the results to the hospital allow the 
doctor at the hospital to give feedback to the personnel in the ambulance so that the 
correct treatment of the patient may start. At the same time, the hospital can prepare 

25 the correct activity for the incoming patient. The system could be very important for 
ECG analyses for all non-specialists in the field if they have to analyse an ECG curva- 
ture. 

The system can analyse drug influence on a number of persons, where analyses are 
30 made before and repeated after drug influence, where selected parameters are com- 
pared and/or combined. It is, hereby, achieved that a drug might be tested for having 
influence on the ECG curvature of a number of persons. This can be very important 
for acceptance of new drugs. This svstem and also described as a method is able to 
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relatively short period, and where the decision if a new drug should be rejected be- 
cause of having negative influence of the ECG, or the drug can be accepted. This deci- 
sion can be taken relatively fast. 

The scope of the invention can also be fulfilled with a method for analysing the drug 
5 influence on the ECG curvature if the method incorporates the steps of: 

receiving ECG curvature from a source, 

- indicating a number of different parameters contained in the received ECG 
curvature, 

10 - storing the parameters in storage means, 

selecting disease specific parameters in the storage means 

- selecting parameters from at least three groups, which groups comprise pa- 
rameters of symmetry, flatness, duration and/or complexity 

- combining selected parameters in mathematical analysing means 

15 - representing the result of the mathematical analysis as a point in a coordi- 

nate system comprising at least one axis, 

- comparing the actual placement in the coordinate system with a number of 
reference parameters stored in a memory, 

- analysing the QT curvature of the ECG for indicating drug induced 
20 changes. 

In this way as already described, a very effective analysis of the ECG curvature is 
achieved. 

25 The analysing process can be repeated in the system for further selected parameters in 
order to achieve more reliable results. Hereby, it is achieved that the system or the 
method can be repeated several times with different combinations of parameters. With 
the system, a deviation of parameters from the stored data indicating symptoms of 
Long QT Syndrome or drug influence may also be interpreted for further reference. 

30 

The system or method analyses the parameters chosen from at least three main groups, 
such as groups of parameters of symmetry, flatness, complexity and duration relating 
to the actual FCG curvature. In this wav, it is achieved that the parameters are grouped 
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30 



cific number of possible parameters. Keeping the number of parameters relatively 
small, the analysis takes place in a faster way. 



The group of symmetry might comprise at least the following parameters: 



SI Symmetry evaluated from Tstart to Tend, calculated by the formula: 

I 

( Tend \$ 

Sl= SCn-m^-wM , 

\n=Tstart J 

where w[n] = v\n\l m 0 , 

Tend 

10 m l = y^W'wfn], 

n=Tstart 
Tend 

m 0 = 

n— Tstart 

and v[n] is the ECG signal. 



iena 



15 S2 Symmetry with Tpeak as mean evaluated from Tstart to Tend, calculated by 
the formula: 

l 

( Tend v 

£2 = — Tpeak) 3 - w[n] 9 

\n=Tstart J 

where w[n] = v[n]/ m 0 , 

Tend 

fn o= 2>M 

n-Tstart 

20 and v[n] is the ECG signal. 



S3 Symmetry with Tpeak as mean evaluated in a symmetric interval of 10% of 

the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 

l 

( Tpeak+0.05{Tend-Tstart) \ 3 



25 S3 = 



y (n — Tpeak) 3 • w[n] 

\ti=Tpeak-0.05(Tend-Tstart) J 

where w[n] = v[n]/ m 0 , 



Tend 



n=Tstart 

and v[n] is the ECG signal. 



S4 Symmetry with Tpeak as mean evaluated in a symmetric interval of 20% of 
the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 
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1 

f Tpeak+O.HTend-Tstart) \3 

S4= ^(n- Tpeak) 3 ■ w[n] 

^n=Tpeak-OMTend-Tstart) J 

where w[n] = v[n]/ m Q , 



Tew/ 



n=Tstart 

and v[n] is the ECG signal. 



S5 Ratio of the time interval "Tstart to Tpeak" and the time interval "Tpeak to 
Tend", calculated by the formula: 
Tpeak - Tstart 



S5 = 



Tend — Tpeak 



S6 Ratio of the average slope from Tstart to Tpeak and from Tpeak to Tend 
S6 = 



Slope TstartTpeak 



15 Slope 



Slope TpeakTeitd 
where 

v[Tpeak] - v[Tstart] 



TstarUTpeak _ 

vjTend]-vjTpeak] 
Mope TpeakJend - Tend _ Tpeak 

and v[n] is the ECG signal. 



20 



25 



S7 Variation evaluated from Tstart to 

l 

( Tend \ % 

S7= £(n-m 1 ) 2 -w[n] , 

\n=Tstart J 

where w[n] = v[rc]/ m 0 , 

Tend 

n-Tstart 
Tend 

n=Tstart 

and v[n] is the ECG signal. 



Tend, calculated by the formula: 



30 SS = 



S8 Variation with Tpeak as mean evaluated from Tstart to Tend, calculated by 
the formula: 

l 

f Tend 

J^(n - Tpeak) 2 -w[n] 



\n=Tstart 

where w\n] = v\n\/ m n , 
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and v[n] is the ECG signal. 
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S9 Variation with Tpeak as mean evaluated in a symmetric interval of 10% of 
the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 



^ Tpeak+Q.05iTend-Tstart) 

S9= ^ (n - Tpeak) 2 ■ w[w] 

K n=Tpeak-0.05'(Tend-Tstart) J 

where w[n] = v[rc]/ ra 0 , 



l 



n~Tsiart 

and v[n] is the ECG signal. 



15 



20 



25 



30 



S 1 0 Variation with Tpeak as mean evaluated in a symmetric interval of 20% of 
the Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 



( Tpeak+0. l{Tend-Tstart) 



510 = 



\n=Tpeak~0. \-(Tend*-Tstart) 

where w[n] = v[n]/ m 0 , 

Tend 

tr h>= S v H 

n-Tstart 

and v[n] is the ECG signal. 



^ (n~ Tpeak) 2 • 



i 



Sll The Hill parameter, K^, evaluated by least square fitting of the repolarisation 
integral, RI(t), from the Jpoint to the following Ponset as described by Kanters et al., 
"T wave morphology analysis distinguishes between KvLQTl and HERG mutations in 
long QT syndrome", Heart Rhythm (2004) 3, 285-292: 

r .» a 

RI(t) = V n 



max 



K 11 -4- f" 



J 



S12 The Hill parameter, K^, evaluated by least square fitting of the repolarisation 
integral, RI(t), from Tstart to Tend analogous to the method described by Kanters et 
al., "T wave morphology analysis distinguishes between KvLQTl and HERG muta- 
tions in long QT syndrome", Heart Rhythm (2004) 3, 285-292: 



f 



RI(t) = V x 



max 



t" 



K K ,n" +t " 



J 



Fl 



Flatness evaluated from Tstart to Tend, calculated by the formula: 

l 

( Tend ^4 

Ml 



WO 2005/058156 PCT/DK2004/000722 

8 



10 



20 



25 



where w[«]= v[n]/ m 0 , 



Tend 



n=Tslart 
Tend 



n-Tstart 

and v[n] is the ECG signal. 



F2 Flatness parameter, F 1 , normalized by the size of the R wave, calculated by 
the formula: 

F2= F1 



v[Rpeak\— v[jpo int] | ' 



where v[n] is the ECG signal. 



F3 Flatness with Tpeak as mean evaluated from Tstart to Tend, calculated by the 
1 5 formula: 



F3 = 



2 

f Tend - -^4 



(n — Tpeak) 4 • w[n\ 

\n~Tstart J 

where w[n] = v[n]/ m 0 , 



Tend 



m<> = Z V W 

n=Tslart 

and v[n] is the ECG signal. 



F4 Flatness parameter, F3, normalized by the size of the R wave, calculated by 
the formula: 

F3 

FA = 



v[i?pea&] - v[jpo int] | ? 



where v[n] is the ECG signal. 



F5 Flatness with Tpeak as mean evaluated in a symmetric interval of 1 0% of the 
3 0 Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 

I 

f Tpeak+Q.05iTend-Tstart) \ 4 

F5 = 5\ ( n ~ Tpeak) 4 • w[n\ 9 

^n^Tpeak-Q.OSiTend-Tstart) J 

where w[n\ ~ v[n]/ m Q , 

Tend 
n-Tstart 

and v[n] is the ECG signal. 

35 
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F6 Flatness parameter, F5, normalized by the size of the R wave, calculated by 
the formula: 

F6= FS 



v[Rpeak] - v[jpo int] | ' 



where v[n] is the ECG signal. 



F7 Flatness with Tpeak as mean evaluated in a symmetric interval of 20% of the 
Tstart-Tend-interval surrounding Tpeak, calculated by the formula: 

I 

^ Tpeak+Q.liTend-Tstart) ^4 

F7= ^(n-Tpeak) 4 -w[n] , 

K n=Tpeak-0. \-{Tend-Tstart) j 

where w\n\ - v\n\l m 0 , 

Tend 

= 2 v[n] 

n = Tstart 

and v[n] is the ECG signal. 



m 0 



F8 Flatness parameter, F7, normalized by the size of the R wave, calculated by 
the formula: 

Fl 



I v[i2pea&] - v[jpo int] | 



where v[n] is the ECG signal. 



F9 Ratio of the total area under the T-wave from Tstart to Tpeak and the corre- 
sponding time interval, calculated by the formula: 

Tpeak 

JTQ — - n— Tstart 



Tpeak - Tstart ' 
where v[n] is the ECG signal. 



F10 Flatness parameter, F9, normalized by the size of the R wave, calculated by 
the formula: 

F9 

F10- 



v[/?peaAr] - v[jpo int] | ' 



where v[n] is the ECG signal. 



Fl 1 Ratio of the total area under the T-wave from Tpeak to Tend and the corre- 
sponding time interval, calculated by the formula: 
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10 
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35 



2 v W 

J n=Tpeak 

Tend — Tpeak 9 
where v[n] is the ECG signal. 



F12 Flatness parameter, Fll, normalized by the size of the R wave, calculated by 
the formula: 

Fll 

F12 = 



v[Rpeak] — v[jpo int] | ' 



where v[n] is the ECG signal. 



F13 Ratio of the total area under the T-wave from Tstart to Tend and the corre- 
sponding time interval, calculated by the formula: 

Tend 

p 1 ^ 3 n— Tstart 

Tend — Tstart 
1 5 where v[n] is the ECG signal. 



20 F\4 = 



F14 Flatness parameter, F13, normalized by the size of the R wave, calculated by 
the formula: 

F13 



v[i?peaA:]— v[jpoint] | ? 



where v[n] is the ECG signal. 



25 F15 Ratio of the T wave height and the T wave width, calculated by the formula: 

_ v[Tpeak] 



F15 = 

Tend — Tstart 



where v[n] is the ECG signal. 



F16 The T wave height, calculated by the formula: 

F16~ v[Tpeak], 

where v[n] is the ECG signal. 



F 1 7 Average slope from Tstart to Tpeak, calculated by the formula: 
F17 = v l- T P eak l- v i Tstart ] 

Tt 7_ nn„ + 5 
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F 1 8 Average slope from Tpeak to Tend, calculated by the formula: 



F18 = 



v[Tend] — v[Tpeak] 



Tend — Tpeak 
where v[n] is the ECG signal. 



10 



F19 The Hill parameter, n, evaluated by least square fitting of the repolarisation 
integral, RI(t), from the Jpoint to the following Ponset as described by Kanters et al., 
"T wave morphology analysis distinguishes between KvLQTl and HERG mutations in 
long QT syndrome", Heart Rhythm (2004) 3, 285-292: 

f » A 



RI(t) = K 



max 



t" 



kKJ +t" j 



15 



20 



F20 The Hill parameter, n, evaluated by least square fitting of the repolarisation 
integral, RI(t), from Tstart to Tend analogous to the method described by Kanters et 
al., "T wave morphology analysis distinguishes between KvLQTl and HERG muta- 
tions in long QT syndrome", Heart Rhythm (2004) 3, 285-292: 



r 



RI(t) = V 



max 



n 



v it 



J 



25 



F21 The Hill parameter, V max , evaluated by least square fitting of the repolarisa- 
tion integral, Rl(t), from the Jpoint to the following Ponset as described by Kanters et 
al., "T wave morphology analysis distinguishes between KvLQTl and HERG muta- 
tions in long QT syndrome", Heart Rhythm (2004) 3, 285-292: 

RI(t) = K 



max 



K K m n +t" j 



30 



35 



F22 The Hill parameter, V max , evaluated by least square fitting of the repolarisa- 
tion integral, RI(t), from Tstart to Tend analogous to the method described by Kanters 
et al., "T wave morphology analysis distinguishes between KvLQTl and HERG muta- 
tions in long QT syndrome", Heart Rhythm (2004) 3, 285-292: 



r 



RI(t) = K 



max 



t" 



J 



QTc The Q-T interval normalized by the square root of the R-R interval according 
40 to "Ra7:ett's formula: 



WO 2005/058156 



12 



PCT/DK2004/000722 



10 



15 



20 



30 



02c — Tend - Qonset 



D2 The time interval from Tstart to Tend, calculated by the formula: 
D2 = Tend - Tstart 



D3 The time interval from Tstart to Tpeak, calculated by the formula: 
D3 = Tpeak - Tstart 



D4 The time interval from Tpeak to Tend, calculated by the formula: 
D4 = Tend - Tpeak 

The group of complexity might contain at least the following parameters: 

CI : Number of local maxima between Tstart and Tend; the minimum number is one. 



C2: Number of phases between Tstart and Tend, where a phase is defined as a singly 
connected part of the wave that is entirely above or entirely below the iso-electric line; 
the minimum number is one. 



25 The parameters previously described can also be calculated and stored as intra- and 
inter-lead means and standard deviations. 



The groups of parameters could contain further parameters, and the groups may con- 
tain a number of subgroups. 



When combining parameters from different groups, a much better result is achieved 
than when only using parameters from the same group. The parameters can be an ele- 
vation of the curve; they can be the morphology of the curve; or they could be time- 
deviations as an example of possible parameters. When combining parameters, a pre- 
35 cise analysis can take place because a specific combination of parameters can indicate 
Long QT Syndrome or drug influence on ECG curvature and it is possible to effec- 
tively select between ECG-signals that look very much alike, but which indicate dif- 
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A selection of these parameters is possible so that special genetic combinations are 
known with reference to the different stored parameters. The system can be updated 
by new data selected from different sources. 

5 The system and/or method can analyse the QT curvature of the ECG for indicating 
Long QT syndrome. This way, the Long QT syndrome can be indicated in an objec- 
tive and effective manner which might occur in postsyncopal cardiac examination. 

The method can differentiate between different genotypes of the Long QT Syndrome, 
1 0 which is important for the treatment. It can, hereby, be achieved that the correct medi- 
cal treatments can be started. The system and the method can be used for test of drug 
influence on ECG curvature. 

The system can be trained, where the parameters' values are calculated for individual 
1 5 subjects, where an analysis of the parameters is performed such as a pattern classifi- 
cation method based on supervised learning, such as Discriminant Analysis, Nearest 
Neighbor Techniques, Multilayer Neural Networks, Decision Trees and Rule Based 
Methods or combinations of these. 

20 The final classification function is at least based on data from at least one LQT or drug 
influenced group and Normal subjects stored as a training set with the consequences 
that the classification method is improved by adding new subjects to the training set, 
which new subject can be tailored to demographic or gender differences. In addition it 
is achieved that reference values based on the training set can be selected from the 

25 most critical group of persons with reference to the parameters that are going to be 
tested. 

Once the parameters' values are calculated for individual subjects the mathematical 
analysis chooses the optimal (small) parameter set out of the complete set (large) from 
30 all categories, which values are stored as ref. values. It should be made clear that the 
final classification functions are based on data from at least one LQT or drug influ- 
enced group and Normal subjects (the training set) with the consequences that the dis- 



WO 2005/058156 PCT/DK2004/000722 

14 

ing set, but also that the method can be tailored to demographic differences (for exam- 
ple in California LQT2 patients might be somewhat different than in Denmark: the 
method can cope with this simply by training the system with people from California 
for use in California and with people from Denmark for use in Denmark) or to other 
5 differences (for example gender differences). Other examples could be age differ- 
ences, difference between infants and adults. 

This invention also comprises the use of a system for analysing ECG curvature for test 
of drugs, which system has input means connected to an ECG source, wherein at least 

10 one among a number of different parameters is isolated and stored in the system, 
where the different parameters of a received ECG curvature are indicated and/or iso- 
lated for indicating possible symptoms, where a number of selected parameters, are 
combined in at least a first mathematical analysis, where the result of the analysis is 
represented as a point in at least one coordinate system, comprising at least one axis, 

1 5 where the system compares the actual placement in the coordinate system with a num- 
ber of reference parameters stored in the system, for indicating symptoms having in- 
fluence on the ECG curvature, and analysing the QT curvature of the ECG for indicat- 
ing drug induced changes to the ECG curvature, where the parameters of the ECG 
curvature are calculated before and after a drug test for a number of subjects, where 

20 the difference for selected parameters between before and after testing is calculated for 
each subject, where a mathematical analysis of selected parameters for a number of 
subjects gives statistical significance for at least one of the following decisions: 
"accept of the drug", 
"rejection of the drug". 

25 "further testing of the drug". 

A very effective way of accepting or rejecting a drug is achieved. 

Below are described one possible method and a system to illustrate the invention. 
30 Abstract 

The Long QT Syndrome is a genetic disorder characterized by abnormal cardiac repo- 
larisation resulting in prolonged QT duration, syncopal episodes and increased risk of 
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than 90% of all LQTS patients. The QT interval duration is the only ECG-based quan- 
tifier of LQTS used in clinical practice today. However duration is only a gross esti- 
mate of repolarisation and does not allow perfect discrimination between KvLQTl, 
HERG and normal subjects. Studies have shown that T-wave morphology parameters 
5 are useful discriminators in LQTS, but no single parameter has proven to be sufficient. 
In this study we present a novel multivariate discrimination method based on a combi- 
nation of T-wave symmetry-, flatness- and duration parameters. 16 subjects were in- 
cluded in the study - 8 normal, 5 HERG and 3 KvLQTl patients. Genotypes were 
known for all LQTS patients, but one. Standard 12 - lead ECG's were recorded on 

10 each subject. An automatic EGG event detection algorithm was implemented. The 
signal was highpass filtered and normalized with respect to the isoelectric level to en- 
sure a stable baseline. 4 parameters describing the duration of repolarisation, 6 sym- 
metry- and 15 flatness parameters were calculated to characterize each of the T- waves. 
The mean values of lead V5 and the interlead standard deviations were used as pa- 

1 5 rameter values. Stepwise discriminant analysis was performed to obtain two discrimi- 
nant functions based on the five strongest discriminatory parameters. The resulting 
discriminant functions include 2 duration-, 2 symmetry- and 1 flatness parameter. The 
two functions classify all subjects correctly (p> 0.0001, p<0.005). Further discrimi- 
nant analysis with a reduced number of parameter categories implied that superior 

20 classification is obtained when using all three parameter categories presented. A com- 
bination of parameters from the three categories symmetry, flatness and duration of 
repolarisation was sufficient to correctly classify ECG recordings from the KvLQTl, 
HERG and normal subjects in this study. This multivariate approach may prove to be 
a powerful clinical tool. 

25 

1. Introduction 

The Long QT Syndrome (LQTS) represents a hereditary genetic disorder characterized 
by the presence of prolonged QT duration on the ECG, syncopal episodes due to po- 
lymorphic ventricular tachycardia (torsade de pointes), and arrythmogenic sudden 
3 0 cardiac death. 

Mutations involving 6 different genes have been identified in LQTS subjects. These 
mutations result in structural and functional changes in ion-channel proteins and cur- 
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repolarization patterns. The most prevalent genes affected in LQTS patients are 
KvLQTl and HERG which account for more than 90% of LQTS genotype patients. 
The current study focuses on carriers of these two genes. Although some attempts 
have been made to develop quantitative measures that link different repolarisation 
5 abnormalities to specific LQTS related channel-opathities these methods have so far 
failed to provide a solid diagnostic yield. In current practice the duration of the QT 
interval is the only widely accepted quantifier of ventricular repolarisation. Yet, it has 
been recognized that the duration of the QT interval is only a gross estimate of repo- 
larisation since T-wave morphology is also important when characterizing the QT 
10 interval. This is evidenced by the fact that approximately 10% of all mutation carriers 
have a normal Bazett corrected QTc (<440ms) and 40% of KvLQTl and HERG carri- 
ers show QTc values between 410-470 ms that overlap with non-carriers. Conversely 
only 2% of all carriers present with a normal ST-T pattern and a normal QT interval. 
Morphological aberrations thus carry major implications for the identification of ab- 

4 

1 5 normal repolarisation and have been included as diagnostic criteria equivalent to that 
of a positive family history for LQTS. 

Studies have shown that affected KvLQTl patients generally show broad based T- 
waves with a normal to relatively high amplitude and often without a distinct T - wave 
20 onset. For individuals with mutations involving the HERG gene the aforementioned 
studies have generally found low amplitude T-waves with bifid T-waves in 60% or 
more of the carriers. 

Cardiologists already include a qualitative assessment of T-wave morphology from the 
25 ECG in order to obtain information that augments the clinically established QT inter- 
val measurement and facilitates discrimination between LQTS genotypes. However 
qualitative description of repolarisation morphology may be biased due to intra- and 
interpersonal variability thus indicating the need for a standardized quantitative meas- 
ure of this parameter. 

30 

In the following is presented a novel multivariate categorization method that allows 
discrimination between KvLQTl, HERG and normal individuals based on Twave 
mornholoffv recorded from 12-lead ECG's. Hallmark morphological features of T- 
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ing three primary T-wave characteristics to be assessed. These characteristics are 
symmetry, flatness and duration. 

2. Methods 

2.1 Subjects 

5 The study included ECG recordings from 8 female and 8 male subjects. The subjects 
were divided into four groups; 3 KvLQTl (aged 20-48, 2 females), 5 HERG (aged 13- 
76, 2 females), 8 normal (aged 23-31, 4 females). Genotypes were known for all 
KvLQTl and HERG subjects with a single exception: 1 patient was categorized as a 
KvLQTl subject by anamnesis and ECG-analysis. In the normal group there were no 
1 0 reports of prior cardiac diseases or LQTS family precedent. 

2.2 Data collection 

Data acquisition was carried out with the subjects resting in supine position. The 
equipment used for data acquisition was a portable digital ECG recording system, 
"Cardio Perfect Resting ECG system" manufactured by Cardiocontrol. Recording was 
15 divided into three sessions. Data was collected from 8 leads (Mil, V2-V6) with a 
sampling rate of 1200Hz. Signal recording length was 75 s. in the first session and 150 
s. in the last two sessions. 

Following data acquisition, SCP files generated by the Cardio perfect software were 
20 exported from a MSDE/SQL7 server and subsequently converted to .MAT files using 
SCP-Batch Converter. 

2.3 Algorithm for detection of events in the ECG 

To facilitate evaluation of the repolarisation process and the QT interval, several 
25 events in the ECG were detected (Qstart, Rtop, Tstart, Ttop and Tend). An algorithm 
for detecting these events was implemented in Matlab 6.0. 

The method is based on prior work published by Laguna et al. and uses adaptive 
thresholding techniaues annlied to a dieitallv filtered and differentiated sienal. A mi- 
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nor extension to the algorithm was incorporated to enable the detection of Tstart 
Tstart was detected with a technique equivalent to the technique for detecting Tend. 
Figure 2 shows an example of the result of the event detection algorithm. 

Figure 2. Important events that are used to describe repolarisation are marked by dots 
5 by the event detection algorithm. The algorithm is able to detect the events on all 8 
recorded leads. 

2. 4 Preliminary signal processing 

Evaluation of the QT interval and the repolarisation process was done on the basis of 
an ECG signal with stabilized baseline. This was achieved through preliminary signal 

10 processing. The "raw" ECG was filtered by a Kaiser window high pass filter with a 
cut-off frequency of 0,5 Hz, 40 dB damping in 0,25 Hz and 0,1 dB ripple in the pass- 
band. Other filters are subsequently used: a lowpassfilter for noise reduction and a 
notch filter for reduction of 50 Hz or 60 Hz interference. The isoelectric line is defined 
as the straight line that connects the PQ interval before the QT interval at hand and the 

1 5 PQ interval after the QT interval at hand. The iso-electric line relative to zero is sub- 
tracted from the QT interval analysed. After filtering, the signal had an almost stable 
baseline. In order to improve stability, isoelectric lines in the signal were estimated 
from one P-Q interval (Qstart minus 20 ms) to the following P-Q interval (Qstart mi- 
nus 20 ms). The signal was then normalized by subtracting the line value from the 

20 corresponding signal values. This process is shown in figure 2. 

2. 5 T-wave morphology parameters 

In order to characterize the T-wave morphology, a number of parameters were se- 
lected. The parameters were chosen to cover each of the three categories: Twave 
symmetry, T-wave flatness and duration. The parameters are listed and described in 
25 table 1. 

Parameters S1-S4 and F1-F8 is based on the calculation of modified skewness and 
kurtosis measures defined as symmetry and flatness in the following. Inspired by the 
summary measures of probability distributions used in the field of statistics the T- 
30 waves were modelled as probability mass distributions (figure 3) and assigned a centre 
(mean), width (standard deviation! an asvmmetrv measure and a convexitv measure. 
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Asymmetry and convexity calculations were then carried out based on the modified 
skewness and kurtosis measures (3 rd and 4 th order moments) as follows: 

The total area under the signal, mO, was calculated: 

N~l 

The signal was normalized by the value of the area, mO: 
w{n] = v[n] I m Q 

1 0 Normalization facilitated the calculation of the moment functions, since w[n] shares a 
fundamental property with the probability mass function: A total area of 1 . 
The 1 st order moment, ml, was calculated, ml is the mean of the signal: 

AM 

H=0 

15 

The 2 nd order moment, m2, was calculated. m2 is the standard deviation of the signal: 

20 Figure 2 Isoelectric lines (dashed lines) in the signal are calculated from one P-Q in- 
terval to the following P-Q interval (Qstart - 20 ms). The line values are subtracted 
from the corresponding ECG signal values giving the distances v(n). The result of this 
procedure is shown as an area plot with basis on the zero-line. 

Figure 3, a) Example probability mass distribution used when calculating standard 
25 skewness and kurtosis measures, b) Modified frequency distribution used in this study 
for calculating the modified skewness and kurtosis measures. Signal values v(n) . are 
shown in figure 2. 
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PARAMETER DESCRIPTION 



Symmetry 

5 1 Skewness evaluated from Tstart to Tend. 

52 Skewness evaluated from Tstart to Tend with Ttop as mean. 

5 S3 Skewness evaluated in a symmetric interval, 10 % of the Tstart-Tend interval 

surrounding Ttop with Ttop as mean. 

54 Skewness evaluated in a symmetric interval, 20 % of the Tstart-Tend interval 
surrounding Ttop wit Ttop as mean. 

55 Ratio of the time interval from Tstart to Ttop and the time interval from Ttop to 
1 0 Tend. 

56 Ratio of the average slope from Tstart to Ttop and from Ttop to Tend. 



Flatness 

Fl Kurtosis evaluated from Tstart to Tend. 
F2 F 1 normalized by the absolute Rtop-Qnadir value. 
15 F3 Kurtosis evaluated from Tstart to Tend with Ttop as mean. 
F4 F3 normalized by absolute Rtop-Qnadir value. 

F5 Kurtosis evaluated in a symmetric interval, 10 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean. 

F6 F5 normalized by absolute Rtop-Qnadir value. 

20 F7 Kurtosis evaluated in a symmetric interval, 20 % of the Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean. 

F8 Kurtosis normalized by the value of Rtop with Ttop as mean. 

F9 Ratio of the total area under the T-wave from Tstart to Ttop and the correspond- 
ing time interval. 

25 F 1 0 F9 normalized by absolute Rtop-Qnadir value. 

Fl 1 Ratio of the total area under the T-wave from Ttop to Tend and the correspond- 
ing time interval. 
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F13 Ratio of the total area under the T-wave from Tstart to Tend and the correspond- 
ing time interval. 
F 1 4 F 1 3 normalized by absolute Rtop-Qnadir value. 
F 1 5 Ratio of the height of Rtop and the width of the Tstart-Tend interval. 

Duration 

QTc The Q-T interval normalized by the square root of the R-R interval according to 
Bazett's formula. 

D2 Time interval from Tstart to Tend. 
D3 Time interval from Tstart to Ttop. 
D4 Time interval from Ttop to Tend. 

The table above shows a Complete list of the parameters used to characterize T-wave 
morphology. Parameters belong to one of three categories: symmetry, flatness and 
duration. 

The 3 rd order moment, m3 5 was calculated. m3 is the modified skewness of the signal: 

^ = (E("- m i) 3 * M W 1/3 

»=0 

Finally the 4 th order moment, m4 ? was calculated. m4 is the modified kurtosis of the 
signal: 



^4=(S(«-^) 4 *^]) I/4 
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2.6 Data analysis in Matlab 

The T-wave morphology parameters for the acquired, pre-processed ECG recordings 
were evaluated using Matlab 6.0. Only valid data were analyzed - i.e. data from leads 
where the signal was not corrupted by high frequency noise and where the event de- 
5 tection algorithm was successful in detecting the relevant events with satisfactory pre- 
cision. Parameter means and standard deviations were calculated for every T-wave in 
the signal on all leads. A great interlead variation in T-wave morphology may be an 
indicator of LQTS. Interlead variance was therefore examined by calculating the stan- 
dard deviation of the lead means for each parameter. 

10 

Only the parameter means from lead V5 and interlead standard deviations were used 
as final parameter values. Hence, for every parameter in table 1 , two parameters were 
calculated - one with index "meanV5" and one with index "std" e.g. FlmeanV5 and 
Flstd. 

15 2.7 Statistical analysis 

In order to characterize and classify data from the three groups (KvLQTl, HERG and 
normal), the evaluated parameter values were processed using discriminant analysis. 
The analysis was carried out in SPSS version 11.5. The objective of the discriminant 
analysis was twofold: finding parameters that most efficiently discriminate between 
20 the groups and reducing the number of variables. Therefore a stepwise procedure was 
used with the Mahalanobis D2 as the most appropriate distance measure. 

The entry/removal-criteria were adjusted in order to reduce the number of variables in 
the discriminant functions to achieve a 1 :3 ratio between the number of variables and 
25 the population size (N=16). The criteria were empirically chosen to be penny = 0.045 
and premoval = 0.09 providing the desired 5 variables in the discriminant functions. 

Figure 4. Scatterplot showing classification of individuals by genotype. Separation of 
groups was carried out by 2 discriminant functions with 5 variables that characterize 
30 repolarisation by computation of symmetry, flatness and duration. 
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3. Results 

The discriminant functions were based on data from all KvLQTl, HERG and normal 
subjects. The 5 parameters included in both discriminant functions are listed in table 2. 

The discriminative efficiency of both generated functions was statistically significant 
5 after inclusion of all 5 parameters (function 1: pO.0001, function 2: p<0.005). 



Variables Entered 



Step 


Entered 


1 


Fllstd 


2 


QTcmeanV5 


3 


S5meanV5 


4 


D4std 


5 


S4meanV5 



1 0 Table 2. Variables used by the two discriminating functions. Stepwise introduction of 
more variables improved the ability of the functions to discriminate between KvLQTl, 
HERG and normal. 

A scatterplot was generated from the discrimination functions and groupings of indi- 
15 vidual genotypes can be seen in figure 4. The dotted lines were read from the SPSS 
generated territorial map and manually added. The lines reflect borderlines where the 
differences between each pair of discrimination functions are zero. All 16 processed 
ECG's were correctly classified and showed at least one discriminatory characteristic 
as defined by the 5 parameters included in the discrimination functions. Cross valida- 
20 tion of both discriminant functions was done with the leave-one-out method and all 16 
subjects were again correctly grouped. Reducing the number of variables resulted in 
misclassified cases due to lack of one or more discriminatory characteristics. In light 

n~ j: ^i^^a +^ KiarfArm 4^it-frV»^v qhqIvcic nf ihp sple.nteH parameters in 
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order to investigate the individual contributions of each variable to the separation of 
the three primary groups of subjects. Extreme values for all parameters were identified 
and the mean was computed. 

5 The result is plotted in figure 5. As expected the extent of interlead flatness variation 
observed in HERG and normal individuals was lower than that found in KvLQTl sub- 
jects. This is evidenced by the Fllstd parameter in figure 5a. When evaluating pa- 
rameter values S4meanV5 and S5meanV5 (figure 5b, d) the extent of asymmetry in 
KvLQTl and normal was generally less than that of HERG individuals. Both 
1 0 S4meanV5 and S5meanV5 are symmetry parameters and asymmetry in HERG indi- 
viduals was augmented in two ways: When bifid T-waves were present the interval 
from Tstart to Ttop was prolonged due to the definition of Ttop used in this study (the 
last highest point on the T-wave). Also, when the initial portion before Ttop was pro- 
longed in HERG individuals better discrimination was possible. Both phenomena 
15 were observed in HERG subjects. Generally the Bazett corrected QTc observed in 
HERG and KvLQTl was higher than that of normal individuals (figure 5e). However 
overlap existed between all three groups preventing separation of the groups by QTc. 
Since no single parameter included in the discrimination functions was able to sepa- 
rate KvLQTl, HERG and normal, we proceeded to investigate the classification effi- 
20 ciency provided by the three primary categories represented by the parameters in the 
functions. This was carried out by generating new discrimination functions using pa- 
rameters from one category only while excluding the other two. Then, from the new 
discrimination functions three additional functions were generated, this time allowing 
the inclusion of parameters from combinations of two categories. Scatterplots illustrat- 
25 ing the results of this analysis are shown in figures 6a-f. The first two functions (figure 
6a) included parameters that characterize the symmetrical properties of the Twave. 
83.1% of the 16 subjects were correctly classified. Arrows in figure 6a indicate the 3 
misclassified subjects. A second discriminant analysis was performed using flatness 
parameters. This resulted in 93.8% correctly classified subjects. Only one subject was 
30 not correctly classified as indicated by the arrow on figure 6b. The misclassified case 
was the same HERG subject incorrectly classified using symmetry parameters. The 
discriminatory efficiency of duration parameters was also evaluated. Discrimination 
*rm1v«i« r«sn1tfiH in Q3_R% correctlv classified subiects. One HERG subject was mis- 
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waves similar to those found in KvLQTl. However the duration parameters failed to 
identify this morphological feature, thus reducing classification performance. 

It can be noted that improved classification was obtained using flatness or duration 
5 parameters versus symmetry parameters and it seemed reasonable to investigate if 
further classification improvement could be achieved using a combination of several 
parameter categories. 

Figure 5. a) Fllstd -Interlead standard deviation of the ratio between the total area 
10 under the T-wave from Ttop to Tend and the corresponding time interval, b) 
S4meanV5 - Lead V5 mean modified skewness evaluated in a symmetrical interval 
surrounding Ttop and corresponding to 20% of the interval between Tstart-Tend. c) 
D4std - Interlead standard deviation of the time interval from Ttop to Tend, d) 
S5meanV5 - Lead V5 mean of the ratio between the time interval from Tstart to Ttop 
1 5 and the corresponding time interval from Ttop to Tend, e) Lead V5 mean QTc. 

Figures 6d-f show the results of three separate discriminant analysis using combina- 
tions of parameters from two categories. It can be noted that classification of subjects 
was perfect in all cases, even when repolarisation duration was not considered (figure 
20 6d). 

4. Conclusion and discussion 

The initial discriminant analysis performed in this study resulted in perfect classifica- 
tion of all KvLQTl, HERG and normal subjects. In table 2 it was noted that the dis- 
criminant functions included parameters from all three categories; T-wave symmetry, 
25 T-wave flatness and duration. This is in agreement with the initial hypothesis that a 
combination of repolarisation duration and T-wave morphology characteristics could 
improve discrimination between KvLQTl, HERG and normal. 

To understand why some subjects were misclassified using a reduced set of parameter 
categories (figures 6a-c) the duration parameters and morphological characteristics of 
30 all 16 ECG's were examined. 
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Using only symmetry parameters, 3 subjects were misclassified. However no obvious 
visual characteristics on the three misclassified ECG's could be identified that ex- 
plained the incorrect classifications. The Bazett corrected QTc was 347ms for the 
normal subject, 425ms KvLQTl, 476ms HERG. Although an obviously prolonged 
5 QTc was present in the misclassified HERG subject it was not identified using sym- 
metry parameters alone. 

Discriminant analysis using parameters from the flatness category resulted in only 1 
misclassification. Again no visual characteristics were identified to account for the 
misclassification. Although it was anticipated that the 

10 

Figure 6. a) The result of discriminant analysis using symmetry parameters resulted 
in three misclassified cases (arrows). Visual inspection of the ECG's revealed no ap- 
parent abnormalities to indicate the reason for incorrect misclassification. b) The re- 
sult of discriminant analysis using flatness parameters. One incorrectly classified 

1 5 HERG subject was identified (arrow) even though no obvious visual abnormality indi- 
cated a different genotype, c) Result of discriminant analysis using duration parame- 
ters. This result illustrates the failure of duration parameters to discriminate between 
KvLQTl, HERG and normal (arrow), d-e) Combinations of parameters from two 
categories illustrate the improvement in classification efficiency when compared to 

20 figures 6a-c evaluation of T-wave flatness would be able to discriminate HERG from 
KvLQTl subjects this was not accomplished by using flatness as a single descriptor of 
repolarization. Performing discriminant analysis based on the QTc parameter as the 
only variable resulted in 1 misclassification. This was not unexpected since it is well 
known that a substantial overlap in QTc values can exist between normal and affected 

25 individuals. The lack of unambiguous discrimination between all groups by use of the 
QTc parameter alone emphasizes the hypothesis that additional parameters are needed 
to classify LQTS individuals. By combining parameters from two categories it was 
found that the discriminatory strength was increased. 

30 (figures 6d-f) This was evidenced by the fact that no subjects were misclassified using 
two categories. A particularly interesting finding, was the perfect separation of all sub- 
jects that was obtained using symmetry and flatness parameters with no duration pa- 
^+^o i^^u^a^a Thic i-^cmH i-miVii^c the* rlicpn'minatorv strength inherent in narame- 
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ters from those two categories. In addition it was found that symmetry or flatness pa- 
rameters combined with duration parameters yielded perfect discrimination between 
all groups. Results from the discriminant analysis using one and two categories indi- 
cate that a combination of more parameter categories strengthen the overall discrimi- 
5 natory power of the classification functions. Combining these findings with the results 
from the three category discriminant analysis initially performed, it is reasonable to 
speculate that a substantially improved discrimination between KvLQTl, HERG and 
normal is possible using all three categories of parameters. 

10 In light of the results obtained in this study we propose a new technique for discrimi- 
nating between KvLQTl , HERG and normal subjects. Through multivariate discrimi- 
nant analysis it was found that a combination of two duration parameters and three T- 
wave symmetry-and flatness parameters was sufficient to classify each of the 1 6 study 
subjects into one of the three distinct groups. Although no single parameter had the 

15 necessary discriminatory strength to classify the subjects, the combination of multiple 
parameters in two discrimination functions was statistically significant (function 1: 
pO.OOOl, function 2: p<0.005). The encouraging results of multivariate repolarisation 
analysis found in this study support the use of symmetry-, flatness- and duration pa- 
rameters to classify LQTS patients. 

20 

The use of the proposed multiple parameter categories to classify KvLQTl and HERG 
genotypes may prove to be a powerful clinical tool in the making. 



